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Abstract 

The theory of second gradient fluids (which are able to exert shear stresses also in 
equilibrium conditions) allows us: (i) to describe both the thermodynamical and 
the mechanical behavior of systems in which an interface is present; (ii) to express 
the surface tension and the radius of microscopic bubbles in terms of a functional 
of the chemical potential; (iii) to predict the existence of a (minimal) nucleation 
radius for bubbles. Moreover, the above theory supplies a 3D-continuum model 
which is endowed with sufficient structure to allow, using the procedure developed 
by Dell'Isola & Kosihski, the construction of a 2D-shell-like continuum representing 
a consistent approximate 2D-model for the interface between phases. In this paper 
we use numerical simulations in the framework of second gradient theory to obtain 
explicit relationships for the surface quantities typical of 2D-models. In particular, 
for some of the most general two-parameter equations of state, it is possible to obtain 
the curves describing the relationship between the surface tension, the thickness, the 
surface mass density and the radius of the spherical interfaces between fluid phases 
of the same substance. These results allow us to predict the (minimal) nucleation 
radii for this class of equations of state. 



Key words: Continuum mechanics ; Gas liquid interface ; Particle size ; Bubbles ; 
Surface tension ; Equilibrium ; Theoretical study ; 
PACS: 47.55.db ; 64.70.Fx ; 68.03.Cd ; 68.03.-g 



Preprint submitted to : European Journal of Mechanics, B/Fluids, 15, 4, pp. 545-568 (1996). 



1 Introduction 



In 1949 Tolman, using the thermodynamical approach due to Gibbs, estab- 
lished a relationship between surface tension and radius of bubbles in equilib- 
rium with their liquid phase [39]. It fails when the bubble has a radius close 
to the nucleation value [29]: indeed by following the ideas of Gibbs [21] it 
is impossible to prove the existence of a 'minimal' nucleation radius, i.e. the 
minimal radius of a bubble or droplet necessary for nucleation in the other 
phase. Cahn and Hilliard [3] proposed more sophisticated models in which the 
nucleation phenomena are accounted for. They introduce a 3D-description of 
the interfacial region, which is characterized as the region in which high mass 
density gradients occur, so that a supplementary term is added to the classical 
free energy expression. This term is in general non-local, i.e. it depends on all 
gradients of mass density [22]. It can be shown however that in many cases 
the first gradientPI of mass density is sufficient to give an excellent approxima- 
tion to the exact solution [17]. This approach was developed in the literature 
mainly using the methods of statistical mechanics [8f9"|16ll42j but also in the 
framework of classical field theories p5l2Uf37] by means of the so-called second 
gradient theory. The second gradient fluids (or capillary fluids as named by 
Casal & Gouin [6]) are fluids which are able to exert shear stresses in equi- 
librium conditions and contact forces concentrated on lines [35]. The theory 
also takes into account the mechanical aspects of the phenomena considered 
in the theories of Laplace, Gibbs and Tolman. Indeed it is easy to prove, for 
second gradient fluids, that the pressure field in equilibrium conditions may 
be non-constant. On the contrary it is an essential feature of the classical 
Laplace theory that the pressure field is spatially uniform. Moreover in the 
framework of second gradient theory it is possible to deduce the PDE gov- 
erning the (equilibrium) mass density distribution from the force balance law. 
Dell'Isola et al. [13] show the consequences of the distinction (first noticed by 
Rocard [34J), between thermodynamic pressure (that deriving from free en- 
ergy) and the mechanical pressure (i.e. the trace of the capillary stress tensor). 
Indeed in the quoted papers it is shown that a definition of the surface tension 
and of the radius of a bubble is possible also in the 3D-approach by using the 
second gradient theory and by means of the concept of the 'equivalent Laplace 
bubble'. In this way it is possible to obtain theoretical predictions for minimal 
nucleation radii and dependence between the surface tension and the bubble 
radius; moreover it is possible to use the 'equivalent Laplace bubble' to study 
the stability of the critical nuclei [21]. Unfortunately in the case of spherical 
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bubbles the ODE determining mass density profile cannot be easily studied, so 
that it is necessary to use numerical simulations in order to be able to make a 
theoretical prediction. We explicitly remark that in the present paper a 'shell- 
like' model for the interface is used, following the ideas developed in [TU]. In 
this approach there is no use of 'Gibbs excesses' quantities (for a detailed 
treatment of these concepts see for instance [7]), so our paper substantially 
differs from [17] in which an analogue analysis for drops is developed. The 
main result of this paper is a new relation between surface tension and the 
radius of bubbles which holds also for microscopic radii and which allows the 
determination of the minimal nucleation radius. The paper is organized as fol- 
lows: in Section 2 we collect the main results due to Gibbs- Tolman, Korteweg, 
Dunn-Serrin, Gurtin and Cahn-Hilliard and then we compare the results with 
experiments; in Section 3 we state the theoretical framework of second gradi- 
ent theory and its application to bubbles; in Section 4 we study the equation 
for mass density profiles for some of the most general two-parameter equations 
of state and solve it numerically; in Section 5 comparisons with experiments, 
together with an explicit prediction about 'minimal nucleation radii' of some 
substances, are carried out; finally some conclusions are collected in Section 
6. 



2 Model of interface between fluid phases: the problem of nucle- 
ation of microscopic bubbles 

2.1 Nucleation of microscopic bubbles 

In [TT][T2"] it is shown that Tolman formula [39J is an universal property of 2D- 
Gibbs models for the interface between phases. In fact the Tolman formula is 
a consequence of the Laplace formula and some simple assumptions (whose 
origin will be found in the works of Gibbs [21]) about the thermodynamics of 
2D-continua. These assumptions are basically: (i) the possibility of defining 
surface quantities such as surface mass density, surface free energy, etc.; (ii) the 
validity of standard thermodynamical relationship between these quantities. 
The Tolman formula predicts that for microscopic bubbles or droplets, the 
equilibrium surface tension will decrease exponentially as a function of the 
equilibrium radius. However very soon after its publication in 1949 Tolman's 
formula was criticized by Lamer and Pound [29] who, founding their analysis 
on nucleation data for very small droplets, found no significant decrease of 
surface tension with the radius. The experimental failure of Tolman theory 
has substantially remained unjustified theoretically until now [18,28j. In our 
opinion two facts are relevant to this problem: (i) the failure of any 'classical' 
theory is a consequence of idealizing the interface as a region without thickness; 
surface quantities would be correctly defined stating a 3D-theory of interfaces 
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and then integrating over the interface thickness in order to define surface 
quantities; (ii) the failure of Tolman formula for droplets/bubbles of very small 
radii cannot be surprising because in the Tolman theory no prediction of 
the 'minimal' critical nucleation radius is made, so the Tolman-Gibbs theory 
cannot explain an essential physical phenomenon. Concerning the point (i) 
many authors have tried to define surface quantities as an integral over the 
interface: the first very crude attempt is due to Tolman himself [38]; many 
authors are able to define surface tension in term of mass density profile, or 
its gradient, by following either the ideas stemming from [3] (cf. also [8|ll6f34| ). 
though these results are applicable only for large bubbles/droplets, or making 
use of the Gibbs concept of a dividing surface [17]; a more rigorous approach, 
using a shell-like approach, is followed in [10J. Point (ii) is closely connected 
to point (i): the description of homogeneous nucleation phenomena is possible 
only in a 3D-theory of the interface because we have to 'form' the interface 
with the critical nucleus starting from a single phase. The idealized model 
of the interface as a simple two-dimensional continuum can obviously fail if 
the radius of the microscopic bubble (or droplet) and the interface thickness 
are of the same order of magnitude. It will be shown below that for typical 
nucleation problem a 2D-theory is not suitable as, in that case, the vapor 
bubble surrounded by its liquid phase is constituted mainly or exclusively by 
the 'interfacial region'. As it is, however, still important to attribute a radius 
and an energy to microscopical bubble, we do it by means of an equivalent 
model of Laplace-Gibbs type. In this way it becomes possible to interpret the 
experimental evidence. 

2.2 Cahn-Hilliard-Korteweg capillary fluids 

The Cahn-Hilliard theory was developed just to account for the nucleation 
or more generally to study the behavior of strongly inhomogeneous fluids. In 
this theory a term depending on the gradient of the mass density is added to 
the classical expression for the free energy valid for homogeneous continua. 
Usually, the additional term is assumed to depend only on the gradient of 
density [3"|I3"T] . Van Kampen has obtained the same results using the ideas 
and methods of statistical mechanics applied to a classical van der Waals 
gas [42] . We remark that in general, gradient theories are in many cases an 
excellent approximation of non-local theory, in which the additional term in 
the free energy depends on all the gradients of the mass density and this can 
be derived only by solving an integral equation [T71T42] . On the other hand 
we note that Korteweg, using the ideas of continuum mechanics, proposed 
some constitutive laws for the stress tensor which allowed him to prove the 
Laplace law [27]. However, as proven by Gurtin, classical Cauchy materials 
cannot show behavior of the Korteweg type, because of the second principle of 
thermodynamics [23] . Dunn and Serrin [T4] proposed to change the expression 
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for the power in Cauchy materials (adding the so-called " interstitial work" ) in 
order to avoid the contradiction shown in [25] . Another possibility could be to 
change the fundamental assumptions about contact forces (Cauchy postulate) 
characterizing Cauchy materials. When interpreted in the language of the 
classical continuum mechanics, the above assumptions about free energy lead 
to an expression for stress tensor which differs from that valid for compressible 
fluids as equilibrium shear components appear in it. The form of force balance 
equation valid for second gradient fluids is found in Germain [20]. However 
controversial, the introduction of second order fluids seems to be able to model 
some interesting problems in a simple way, as underlined by Casal, who named 
them "capillary" fluids [I]. We remark that Cahn and Hilliard do not consider 
the mechanical aspect of the nucleation phenomenon, which is accounted for in 
the Laplace approach. For this reason they only consider the thermodynamic 
pressure (i.e. the spherical part of the stress tensor deriving from the classical 
expression of free energy) while Casal shows, for second gradient fluids, the 
existence of a capillary non-spherical stress tensor whose trace includes, but 
does not reduce to, the quoted thermodynamic pressure [5]. Finally we remark 
that Cahn and Hilliard introduced, using their 3D-approach, an expression 
for surface tension valid only for plane interfaces, and never investigated the 
relationship between surface tension and radius (see also [SUS])- 

2.3 Experimental evidence 

From an experimental point of view the failure of Tolman formula was verified 
by Fisher and Israelachvili [18J, though their results deserve some comments: 
their fundamental result is the observation that no significative variation of 
surface tension was observed down to radii near to the molecular diameter. 
However their experimental apparatus cannot reproduce a process of homo- 
geneous nucleation because at such small radii the adhesive effects of their 
'hyperboloidal' bubble can significantly change the boundary value problem 
for the density profile. On the other hand the 'minimal' nucleation radius 
predicted from the Cahn and Hilliard theory or by second gradient theory ap- 
pears typically larger than the molecular diameter by a factor of two or three, 
so reproducing the data typical of nucleation studies (see Section 5 below). 
Another problem connected with the Fisher and Israelachvili experiment is 
that radii are obtained via the Kelvin relation that (though surely 'verified 
down to 4 fim' as quoted by the authors) needs some correction near to the 
'minimal' nucleation radius. We think that it is possible to reinterpret their 
data by using the equivalent model described in the next section. More re- 
cently Kumar et al [28J proposed a slight modification of Tolman results that 
can account for some of the observed nucleation data: a comparison of second 
gradient theory with their values is reported below in Section 5 for some data; 
here we note that the Kumar theory cannot predict a 'minimal' nucleation 
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radius for droplets (the 'critical' radius data in their paper will be interpreted 
more correctly as an unstable 'equilibrium' radius [JT] • The reason is that their 
theory is essentially a modification of the 2D-Laplace-Gibbs theory in which 
the interface has no thickness. 



3 Shell-like interfaces by second gradient theory 



3. 1 Equilibrium equations for a fluid endowed with internal capillarity 



The second gradient theory, conceptually more straightforward than Laplace's 
theory, can be used to construct a theory of capillarity. We aim to study the 
equilibrium properties of a system with an interface under fixed temperature 
(isothermal conditions). In the following whenever stated capital letters P, 
W, etc. denote classical (per unit volume if the quantity is extensive) ther- 
modynamic quantities; small Greek letters e, p denote classical mass-based 
thermodynamic quantities. For example, if ip is the classical mass free energy, 
the volume free energy will be W = pip, and the classical chemical potential 
p is given by: 

dW 

M = -jj- (3-1) 

In the present text the only addition is an internal mass energy e that is a 
function of the density p as well as Vp. The internal mass energy characterizes 
both the compressibility and capillarity properties of the fluid, independently 
of the bodies with which it is in contact. For an isotropic fluid, it is assumed 
that: 

e = e(p,P p ), (3.2) 
where (3 = (Vp) 2 . The equation of equilibrium is written: 

divS-pVfi = 0, (3.3) 

where Q is the extraneous force potential and S is the general stress tensor: 

S = -pI-A(Vp)® (Vpf, (3.4) 

with A = 2pe 5j a p and p = p 2 e iP — pdiv(AVp). Following Rocard, this last 
expression for p is what we call the mechanical pressure [3i 



2 We remark that in [T7] only the first term p 2 e :P , corresponding to the classical 
thermodynamic pressure P, enters in the Laplace law, cf. Eq. (4.9) of the quoted 
reference; however in appendix A the authors agree that the mechanical pressure 
needs consideration. 
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The equation of equilibrium (3.3) is obtained in the clearest manner by using 
classical methods of the variational calculus, i.e. using the virtual work Prin- 
ciple due to d'Alembert-Lagrange (cf. Serrin [HE])- It appears that a single 
term containing A accounts of the second gradient effects in the equilibrium 
equation. The quantity A, like e, depends on p and Vp (see also the discussion 
of the predictions of statistical mechanics on the p dependence in Section 4 
below and references cited therein) . In a study of surface tension based on gas 
kinetic theory [34], Rocard obtained the same expression (3.4) for the stress 
tensor, but with A constant. If A is constant, the internal energy is simply 
written: 

A 

2p 



e(p,/3 p ) = e (p) + —/3 p , (3.5) 



i.e., the second gradient term — (3 p is simply added to the energy eo(p) of 

Zp 

the classical compressible fluid. The thermodynamic pressure for the fluid is 
P = p 2 £o,p, which gives for the mechanical pressure p: 

p = P-\(±P p + p A p y (3.6) 

where A is the Laplace operator. This pressure appears in the boundary con- 
ditions [35]. For P, Rocard uses the van der Waals pressure or two-phase state 
laws. 

In isothermal equilibrium conditions, Eq. (3.3) can also written as: 

p(p) — \Ap + Q = const. , (3.7) 

where p is given by Eq. (3.1), i.e., is the chemical potential, relative to the 
fluid in an homogeneous state with mass density p at a fixed temperature T. 
If £1 is negligible the equilibrium of a liquid bubbla[El surrounded by its liquid 
phase with density p { is represented by a spherically symmetric density profile 
satisfying: 

A (p irr + -p^j = p(p) - p(p[), (3.8) 

We call Eq. (3.8) the Density Profile Equation (DPE). To the DPE we have 
to add boundary conditions: 

(i) because of spherical symmetry the derivative of p vanishes at the origin; 

(ii) as we have assumed that the bubble is surrounded by an homogeneous 
liquid the same condition holds at infinity^* 



3 The same approach could be followed also in the case of vapor droplets, but here 
for sake of brevity we limit ourselves to develop the theory for bubbles. 

4 If pf is the density of the liquid phase in the equilibrium state with plane interface 
and if the constitutive function p is suitably regular, the theory of Fuchs equations 
implies that: 

(iii) for every p l in an interval (p spi , pf), where p spi is the spinodal density value, the 
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Equation (3.8) was formulated by Rocard [33] and by Cahn & Hilliard [3] 
(cf. also Blinowski [2], Truskinovsky [30] and for a more general mathematical 
treatment Peletier & Serrin [3D]). 

With W, the volume free energy being a primitive of p, by multiplying the 
DPE by p r and integrating between zero and r we obtain: 

^p 2 r + 2A f ^rfr = W(p) - W{p v ) - p{ Pl ){p " PJ- (3-10) 
2 ' Jo r 

Moreover by multiplying the DPE by r and integrating again between zero 
and r, we obtain: 

Xrp r (r) + A(p(r) - p J = / [p(p) - n( Pl )]rdr, 

Jo 

and by taking the limit for r — > oo, rp r (r) tends to zero and we obtain: 

rod 

KPi - Pv) = / Hp) - v{pi)Vdr. (3.11) 
Jo 

Such expressions as Eqs. (3.10)-(3.11) can be useful to test the numerical 
calculations of the interface profile, especially the second one because it relates 
the density jump across the interface directly to a simple integral involving 
chemical potential. 

3.2 Nucleation energy of bubbles 

In the theory of Laplace-Gibbs we have for bubbles of radius R and surface 
tension a: 

w = AnR 2 a + ^nR 3 (w(p v ) - W( Pl ) - »( Pl )(p v - Pl )), (3.12) 

with P = pp — W being the thermodynamic pressure; the isothermal equilib- 
rium conditions: 

P(pl) = MPJ; 
P{ Pl ) - P( Pv ) = ~ 

DPE and conditions (i)-(ii) determine uniquely an increasing mass density profile 
p(r) and, in particular, the value p v at the origin; we will say that p{r) satisfies the 
capillary fluid version of Gibbs Phase Rule; 
(iv) p[r) is twice differentiable at the origin; 

v) as (jij^J > 0, p{r) converges at least exponentially to pi when r tends to infinity. 
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are valid only in the context of Laplace theory and transform Eq. (3.12) into: 

4 



w = — nR 2 a. 



(3.13) 



One can conclude that the nucleation energy of the bubble is one third of the 
creation energy of its interface. 

We will now extend this result to the theory of second gradient fluids: in the 
theory of second gradient fluids, by following Cahn and Hilliard [3], we have 
for the nucleation energy of a bubble in a domain D (cf. also Dell'Isola et al 
El): 



w 



D 



W(p) - W{pi) - p( P/ )(p - Pl ) + ^(Vp) 2 



dv. 



(3.14) 



By denoting ip(p) = W{p) - W(p t ) - n{pi){p - p t ), by multiplying Eq. (3.10) 
by r 2 and by integrating the result over [0, +oo] we get: 



A 



r 2 dr + 



V>(pJ + 



2A 



o r 



p 2 r dr 



r dr = 0. 



By integrating by parts and using Eq. (3.10) again, this equation becomes: 



P 2 r+^(P) 



r dr + 



Because of properties of the density profile quoted in Section 3.1, the second 
term vanishes, and we get: 



W = 47T 



Hp) + \p\ 



4 roo 

r 2 dr = -ii \p 2 r 2 dr. (3.15) 
3 Jo 



By denoting R the mean value of r 2 with respect to the measure p 2 r dr, the 
above equation reads: 



w = -txR 2 / Ap 2 dr. (3.16) 
3 Jo 

For a large bubble, i.e., when p z tends to the plane interface value, R represents 

roo 

the radius. The surface tension for a plane interface is / Xp 2 dr l3f8lfT6l . 

Jo 

Therefore Eq. (3.16) extends Eq. (3.13) to microscopic bubbles and reduces 
to it in case of large bubbles. 



3. 3 Comparison of Laplace and second gradient theories. Equivalent Laplace 
bubbles 



In the second gradient theory the stress tensor in the centre of a spherical 
bubble takes the value S = — p v l where p v = P(p v ) — Ap^ Ap^ p=p . The value 
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Pi of the mass density in the liquid phase is attained asymptotically while the 
stress tensor takes the value S = — p/I where pi = P(pi). As the DPE implies 
^P v &p p=Pv = p v [p{p v ) - p(pi)] we have: 



Pv-pi = W(p l )-W(p v ) + p(pi)(p v -Pi). 



(3.17) 



Let us notice that this difference is not equal to the corresponding difference of 
the thermodynamical pressures as, for microscopic bubbles, p(pi) differs from 
p(p v ). As the experimental results (see for instance Fisher & Israelachvili [18J ) 
deal with measurements of stresses, then we have to use p v — pi instead of 
P(p v ) — P (pi) in the comparison between Laplace theory and second gradient 
theory. We can now define the surface tension and the radius of a bubble by 
identifying the nucleation energies and the pressure differences computed in 
both theories as follows: 



= £ and 
which leads to the relations 



4 4 r°° 

-■nR 2 o — -7r / Xp 2 r r 2 dr, 

3 3 Jo 1 



R 



2A / p z r r 2 dr 
Jo 



[W(Pi)-W(p v )-p(pi)(pi-p v )}-^ 



(3.18) 



and 



a 



4 Jo 



p 2 r r 2 dr 



Mpi) -W(p v ) - p{pi){pi- pj]3 



(3.19) 



Equations (3.18)-(3.19) define the radius and the surface tension of a micro- 
scopic 'equivalent Laplace bubble' in our 'shell-like' theory. We remark that 
they differ from Eq. (4.9) of Falls et al [17]; moreover, in Falls et al, R is not 
specified in terms of the density profile. Indeed Figures 18 and 19 of the quoted 
reference are drawn by means of the 'Gibbs excesses' approach, so that these 
plots have in abscissa the Gibbsian 'radius of tension' or 'surface of tension' 



3.4 Thickness and surface mass density of interfaces 



One of the fundamental problems of interfaces is the study of their thickness. 
We note that in the Laplace-Gibbs model the interfaces have no thickness. 
In the above proposed model, in which A is constant, the interface thickness 
is infinite. This is a consequence of the fact that the limit density in the 
phases is reached asymptotically. On the other hand the experiments show 
that the interfaces generally have a very small thickness, i.e., typically a few 
molecular diameters. The thickness takes a macroscopic dimension when the 
critical point is approached or the density is almost uniform (near the so-called 



10 



spinoidal limit [3]). Therefore a sort of coherence length for the interface is 
needed in second gradient theory. Some natural expressions can be proposed to 
characterize the thickness of the interface taking into account only the region 
in which the strongest gradient occurs. The most simple of these expressions 
is the quadratic variance associated to a measure p 2 r r 2 dr (we note that the 
first order variance is zero). We propose: 

/ (r - R) 2 p 2 r r 2 dr 
S 2 = ^ joo : • (3.20) 

2 2 t 

p r r dr 



We call 5 the 'equivalent' thickness of the interface. We remark that there is 
also a natural definition of the bubble radius: 

/ rp 2 r r 2 dr 

Ri = -yoo • (3-21) 

/ p 2 r 2 dr 
Jo 

Numerical simulations show that, in all the cases studied below, Ri and R are 
practically the same (differing by less than 1% except very near to the critical 
point). In the 2D- interface models the notion of surface mass density p a is 
introduced. In the treatment proposed in [10J it is assumed that: 

f-B+S/2 

p a = / pjdr, (3.22) 

JR-S/2 

where j(r) = 1 — 2H(r — R) + K(r — R) 2 is the curvature dependent Jacobian 
pertaining to the change of variables between surfaces parallel to the interface 
in the 2D-models. 

In the next Section we study by numerical simulations, for different two- 
parameter equations of state, the dependence of the interface thickness and of 
the surface mass density as functions of the radius and the temperature. 



4 Qualitative analysis and numerical solutions of the DPE 



4-1 General qualitative analysis 



For the normalized density p(r), the DPE (3.8) reads (cf. [2"f3"] ): 

2 

p rr + -p, = fi(p)-H 0O . (4.1) 

r ' 

Now, subscript means that the derivatives are taken with respect to the nor- 
malized length variable r: in the whole section lengths and densities are nor- 
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malized with respect to C a = \ t — =- and critical density p . The quantity 

V X Pc Z c 

p is the normalized chemical potential and /i M = p{pi)- The explicit normal- 
ized form of the free energy W is given below for some typical equations of 
state. Z c and P c are the critical compressibility ratio and pressure. We will 
look for a solution of the DPE verifying the boundary conditions considered 
in Section 3.1, i.e., p r (0) = p r (oo) = 0, that represents a physically possible 
density profile in the phase transition region describing the equilibrium of a 
bubble (droplet) with its liquid (vapor). The qualitative analysis is based on 
the 'mechanical' interpretation of Eq. (4.1) [12]. This may be regarded as the 
equation of a 'particle' of mass 1 moving in the (//^-dependent) potential: 

U(p,p 00 ) = -(W(p)-p 00 p), (4.2) 

with the 'viscous' time-dependent force (2/r)p r . In the following we refer to 
the 'motion' of the particle as a solution of the Eq. (4.1) starting from r = 0, 
i.e., at the time zero, with zero velocity (cf. Section 3.1) and with a given 
initial density, i.e., a given initial position, p(0). 

It is evident that a fundamental role in our problem is played by the potential 
U(p, /ioo). For the sake of brevity, we limit ourselves to some general consid- 
erations valid for any potential. We assume under very general assumptions 
that: 

(a) for any temperature below the critical temperature there exists a range of 
values for p^, say (p s c pi ' bub , p sp > hdrop ) , for which the potential (4.2) has three 
extrema Pi^p^) < p m (Poo) < P2(poo)> Pi an d P2 are the maxima and p m is 
a minimum; 

(b) for p^ = [i s P i l > bub (or p^ = p s £ l ' drop ), at the so-called spinoidal limit [4"T| . 
the potential has a maximum only in p l (or p 2 ) and an inflexion point in p 2 
(or Pl ); 

(c) the value of the chemical potential at the equilibrium between phases with 
planar interface (saturation) p s ^ lies in the interval (p s c pi,bub , p s c pt ' drop ); for 
Poo = Poo* the values of the potential is the same at the two maxima 
(Maxwell rule), i.e., U( Pl ) = U(p 2 ); 

(d) for p s c pi,bub < p^ < p s ^ the value of the potential in p 2 is always lower than 
that in p x , i.e., U(pi) > U(p 2 ) (bubble case); if p s £ < p^ < p spi ' drop , the 
opposite happens: the value of thepotential in p 2 is higher than that in Pl , 
i.e., U(p 1 ) < U(p 2 ) (droplet case) lj. 

From assumptions (a)-(d) is clear that a 'motion' starting with initial position 

5 The simplest model is obtained by developing the free energy functional, near 
the critical point, as a function of p; the result is W(p) = (p — a\) 2 {p — (I2) 2 (with 
a\ and 02 constants corresponding respectively to the homogeneous gas and liquid 
density), which satisfies the properties (a)-(d) as it can be simply verified. 
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p(0) in the region just at the right (left) of the first maximum p l (second 
maximum p 2 ) of the potential for the case of the bubble (droplet) rolls down 
to the minimum p m and then goes up towards the second maximum p 2 (first 
maximum p x ) where three different cases are possible: 

1) The velocity is insufficient to reach the second maximum p 2 (first maximum 
Pi)\ this happens because of the dissipative term in Eq. (4.1). In this case 
the 'particle' reverses its velocity before reaching the second maximum (first 
maximum) and falls towards the minimum where it will be trapped with 
damped oscillations when r tends to oo; 

2) The velocity is sufficient to reach the second maximum (first maximum); 
then the 'particle' arrives in general with a finite velocity at this point, then 
it falls down in the region of high density p > p 2 (low density p < p-J after 
the maximum, its velocity grows and the density reaches the co-volume 
(zero) limit in a finite time; 

3) The velocity is exactly that one sufficient to reach the second maximum (first 
maximum) (separatrix solution), so that the 'particle' reaches the second 
maximum (first maximum) in an infinite time, i.e., for r tending to oo; at 
the maximum the velocity p is zero, so that the boundary conditions stated 
in Section 3.1 are satisfied. 

We will call pj the value for p at which the potential attains the lower maxi- 
mum, i.e., p f = p 1 X p s £ < p^ < pr ,drop , and p f = p 2 if p^' feufe < p^ < p^*. 
Obviously only the case 3) is physically significant. The force acting on the 
'particle' is zero when r tends to oo as the maximum is an extremum of the 
potential. From Eq. (4.2) it follows: 

KPf) = Poo = MPoo)- 

This implies that (we assume that locally for large [small] values of density 
the function p(p) is invertible), Pf — p^- 

So the solution of the DPE is simply found if we examine the solution of the 
equivalent mechanical problem when starting a 'motion' with an initial value 
p(0) such that we have exactly the separatrix solution. So for a fixed p^ the 
only unknown quantity is p(0) we will find when we solve numerically Eq. 
(4.1). 

We note that far from saturation the equality of (classical) chemical potentials 
is not satisfied because in Eq. (4.1) we do not have zero on the RHS when r 
tends to zero. So for a small bubble an additional contribution due to curvature 
appears in the Gibbs energy if we want to preserve the classical relationship. 

Near the saturation limit the problem becomes substantially different: in the 
case of bubbles the 'particle' starts its motion in the neighborhood of p 1 and 
tends to reach pj = p 2 ; its velocity has to be relatively large and therefore dis- 
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sipation could be very large and inhibit the approach to the second maximum. 
However the system can react in a different way as a full separatrix solution 
linking the two maxima exists exactly at p^ = //^f . In fact at saturation the 
two maxima have the same height (Maxwell rule) and the separatrix starting 
from the first maximum and ending over the second is the correct solution 6 1. 
This is exactly the Van Kampen solution for the planar interface |42j . 



4-2 Free energy and A value in some simple models 



The free energy functions can be chosen for the problem of phase transition 
by modelling them on an appropriate equation of state for the single fluid. 
The properties (a)-(d) will be shared by many approximate models of the 
long range effects such as the van der Waals or the Berthelot equations of 
state. In the following we study only the class of two-parameter equations of 
state though the method is applicable also to more general equations. This 
class includes several models [M] with the standard expression a/V 2 for the 
internal pressure and also some more sophisticated models involving further 
correction of the internal pressure [31] . However the form of the free energy is 
substantially the same leading only to some relevant but quantitative changes. 
We assume that the equation of state can be written in a general form as: 

P(p,T) = ptf(f3p)-ap 2 g(Pp), (4.3) 

where t is the reduced temperature, i.e., t = T/T c ; a = ^ c (with TZm = 

1Z/M., 1Z the gas constant, M. the molecular weight) and (3 = p c b are the 
generalized normalized van der Waals coefficients (which can depend on tem- 
perature T); / and g are suitable independent functions. In Table 1 we list 
the equations of state studied in this paper, the form of functions / and g and 
the values of the normalized coefficients a and (3. The normalized volume free 
energy, derived from the above equation of state, is: 



W(p) = P F (t) -pt f ° f{fip) ^-pa [ P0 g((3p) dp, (4.4) 
Jp p J P 



where p is a reference state at which free energy is p Fo(t). With regards to 
A, Evans [16] shows under very general assumptions that: 

Hp) = -J^- J c[r;p]r 4 dr, (4.5) 



6 This happens because the dissipation, i.e., the term (2/r)p r , cannot play any role 
because the escape time from the first maximum is now infinite and this term tends 
to zero as r tends to infinity. 
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where is the Boltzmann constant and c[r; p] is the direct correlation func- 
tion of the fluid. This implies that A is a functional of the field p. However 
under some degree of approximation [T5] an expression for the direct correla- 
tion function is given by the Percus-Yevik equation (PY-equation) 

V 

c[r;p] = g[r;p)(l - e 

where g[r; p] is the pair correlation function and v is the intermolecular po- 
tential. The expressions c[r] and g[r] can be numerically found in this case 
[15] . In order to evaluate A, we make the drastic approximation of taking 
g[r] = 0(r — ro), where is the Heaviside function and tq is typically of the 
order of magnitude of the molecular radius. 



State Equation 


a 




/ 


9 


van der Waals 


9 
8 


i 

3 


l 

l-x 


1 


Rocard 1 


1.254 


0.430 


(l-f) u 


1 


Rocard 2 


1.319 


0.478 


(l-f) a 


1 


Rocard 3 


1.500 


0.617 


e — X 


1 


Rocard 4 


1.370 


0.515 


i ^ 


1 


Berthelot 


9 1 
8 t 


i 

3 


1 

1-x 


1 


Eberhart-Schnyders 


9 1 

8 1 

42 


1 
3 


1 

1-x 


1 


Peng-Robinson 


1.489(1 - fc(l - t 2 )) 


0.253 


1 

1-x 


1 ! 1-X-V2 
2V2 X-1-V2 



Table 1. Scaled parameters a and (5 and functions / and g for the equations of 
state considered in this paper 
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Developing the Boltzmann factor in the PY-equation to the first order in v, 
this approximation leads to the standard results of Cahn and Hilliard [3] or 
equivalently of Van Kampen [12] (cf. also Gouin [23]) in which: 



We note that the normalizing length C a , can be related to A in the following 
way. In fact we can write C a = y/otC where a = apj RT C and (cf. [3]): 



We note that this last length depends strongly on the intermolecular potential: 



for example for the 12 — 6 Lennard- Jones potential, we have £ = y ll/7r in 
which r = r LJ is the Lennard- Jones radius. 

4-3 Numerical results 

We integrate Eq. (4.1) by a Bulirsh-Stoer integrator [33] starting from an 
arbitrary guess value of p(0) with p r (0) = 0. The values of p(0) are then 
adapted with a sequence of integrations until a good approximation to the 
separatrix is reached. Near saturation, the accuracy of the method is limited 
by the machine precision because the starting density differs from the first 
maximum density p x by progressively smaller quantities. We were limited to 
the usual double precision (10 -16 ) of our computing device. We find that this 
precision is sufficiently high to catch the main features of the approach to 
saturation. To test the procedure we have verified that, starting from the last 
point and going back toward small radii the 'particle' follows almost exactly 
the same motion (backward in time). As a further test we have verified that 
integral relation Eq. (3.11) is satisfied within the error introduced by numerical 
integration. 

Typical density profiles are shown in Figure 1 for the van der Waals equation 
of state. The two parts of the nucleation phenomenon, i.e., the homogeneous 
nucleation and the growth of bubble approaching saturation, are clearly seen 
in this figure. The density profile tends to be diffuse towards the spinoidal 
limit so that the density jump p^ — p(0) goes to zero (cf. curves (a), (b) and 
(c) in Figure 1). Approaching saturation, a region of almost-constant density 
appears before the phase transition region, i.e. a bulk vapor phase, and the 
density jump is practically constant (cf. curves (f), (g) and (h) in Figure 1). In 
Figure 2 we show the dependence of a on R (as defined in Section 3.4) numer- 
ically evaluated at the reduced temperature t = 0.85 for some of the equations 




(4.6) 




(4.7) 
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Fig. 1. Density profiles for van der Waals equations of state for different values 
of //qq (or p^) at t = 0.85 obtained by numerical solution of the DPE. Towards 
saturation, cf. the curves labeled from (e) to (h), the bubble clearly develops a bulk 
vapor phase in its interior. Near the minimal nucleation radius, cf. the curve (d), the 
bulk vapor phase disappears. In the spinoidal limit, cf. the curves from (c) to (a), 
the radius of the bubble again grows but no bulk vapor phase exists and the density 
jump rapidly decreases until it disappears at the so-called spinoidal composition in 
which only the liquid phase exists. 

of state listed in Table 1. The results are qualitatively similar: for small val- 
ues of Poo near the spinoidal value fi b £ b ' spi (i.e. the lower part of the curve 
in Figure 2), a initially vanishingly small, grows slowly whereas the radius 
decreases towards a minimal value that we identify as the 'minimal nucleation 
radius' R m . At R m the surface tension suddenly increases and the density at 
the center of the bubble reaches the range of the vapor phase (cf. also curves 
(d) and (e) in Figure 1). Above the minimal nucleation radius, a approaches 
within a few minimal radii the planar interface value and, as shown by Fisher 
and Israelachivli experiments, remains substantially constant [18]. We note 
that quantitatively, the van der Waals equation of state gives the lowest value 
of planar interface surface tension, the Berthelot and the Peng-Robinson give 
the highest values. The behavior for other values of temperature is qualita- 
tively the same. We note also that the curves end at different radii as they 
approach saturation; these radii have to be considered as equivalent because 
the 'distance' from saturation, measured by the difference between the first 
maximum of the potential and the density at the center of the bubble, is in all 
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Fig. 2. Normalized surface tension a ja n where a n = p c RmT c C vs normalized equi- 
librium radius R/C of microscopic bubbles at t = 0.85 for the following equations of 
state: van der Waals □; Rocard 4 Berthelot x; Peng-Robinson A; Eberhart-Schny- 
ders +. The end-point of each curve toward the saturation have to be considered 
equivalent in the sense of saturation proximity (cf. the text). 

cases less than 10~ 16 . It could be implied by this circumstance that different 
equations of state reach saturation at different radii. Moreover we note that 
curves corresponding to radii between the spinoidal and the minimal critical 
radius values can be considered as lacking physical meaning because, no phys- 
ical bubble can exist in this regime, except fluctuations of typical radius R 
(cf. [3]). So the 'minimal nucleation radius' R m marks the last values, before 
the spinoidal, to which one can give physical meaning. 

In Figure 3 we plot the (equivalent) interface thickness denned in Eq. (3.20) 
as a function of R for some of the different equations of state, always at re- 
duced temperature t = 0.85. At the start of nucleation (spinoidal limit), 6 is 
very large (diffuse interface) and of the same order of magnitude as R: density 
fluctuations have a unique length scale. When R has reached the minimal nu- 
cleation radius, 5 quickly tends to a practically constant value with very small 
deviation from the limit value at saturation. Again all the studied equations 
of state have the same qualitative behavior. Moreover, we note that the curves 
between the spinoidal limit and the minimum radius are essentially the same 
for all the equations of state used here (cf. Figure 3). 

We note that no qualitative differences appear when the temperature is changed. 
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It is also interesting to note that the ratio (5/R) m is, within the numerical 
accuracy, a constant independent of t: its value is reported in Table 2 for all 
equations of state. Moreover this constant is practically the same for all the 
equations of state studied here, i.e., (S/R) m ~ 0.42 — 0.44. The justification of 
this phenomenon can be that, the shell-like quantities as surface tension, ra- 
dius, etc., being defined through an integral over the thickness of the interface, 
the differences between the equations of state will be smeared by integration, 
giving similar results. 



State Equation 


(S/R) m 




p£7p? 


van der Waals 


0.427 


0.771 


0.716 


Rocard 1 


0.443 


0.778 


0.712 


Rocard 2 


0.435 


0.785 


0.739 


Rocard 3 


0.443 


0.768 


0.716 


Rocard 4 


0.432 


0.795 


0.773 


Berthelot 


0.426 


0.768 


0.710 


Eberhart-Schnyders 


0.428 


0.761 


0.683 


Peng-Robinson 


0.439 


0.750 


0.767 



Table 2. Numerically evaluated ratios (5/R) m , Soo/5 m and for the equa- 

tions of state reported in Table 1. The ratios are substantially independent of the 
temperature in the limit of the algorithm error introduced in the numerical solutions 
of the DPE equation. 

In Figure 4 we plot the surface mass density versus R at t = 0.85 obtained by 
using Eq. (3.22) for the same equations of state than in Figures 2-3. We see 
that surface mass is an increasing function of R going towards the spinoidal 
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Fig. 3. Normalized thickness 5/C of the interface vs normalized equilibrium radius 
R/C of microscopic bubbles at t = 0.85 for the following equations of state: van der 
Waals □; Rocard 4 *;Berthelot x; Peng-Robinson A; Eberhart-Schnyders +. 

limit. This behavior is obviously due to the increase of 5 in the same limit. 
Again we note that no qualitative differences appear when the temperature is 
changed. Moreover the ratios 5oo/6 m and are independent of t. Their 

values, listed in Table 2 for all the equations of state, differ by less than 15%. 



5 Minimal nucleation radii: theoretical predictions and comparison 
with experiments 

The comparison of the above results with the experimental values reported 
in literature is a rather difficult task because information for many data can 
be only indirectly recovered. Moreover, as shown in Section 4.2, Eqs (4.5- 
4.7), all the predictions are critically dependent on the parameters of the 
intermolecular potential; in our case a fundamental role is played by £ that 
sets the length scales. Again the correct prediction of surface tension is also 
not simple in the sophisticated microscopic statistical theories (cf. the data 
reported in [16]) because, especially at very small scales, the details of the 
intermolecular potential, and so the particular molecular micro-structure, are 
relevant in the determination of the exact values of surface tension. We prove 



20 



10 



T 



T 




1 




-1# -Q -B 











5 



10 



15 
R/L 



20 



25 



30 



Fig. 4. Normalized surface mass density p a j (p c £) vs normalized equilibrium radius 
R/C of microscopic bubbles at t = 0.85 for the following equations of state: van der 
Waals □; Rocard 4 *; Berthelot x; Peng-Robinson A; Eberhart-Schnyders +. 



in this paper that rather good predictions of surface tension are also possible 
in the framework of the second gradient theory. 

Moreover we are able to produce predictions of nucleation radius. In particular, 
by evaluating R m , we establish a lower bound to the experimental nucleation 
radii reported in the literature. We remark that very few data are found in 
the literature concerning bubbles (cf. [4~I]). so that we are obliged to compare 
these predictions with data for the droplets found in Lamer and Pound or 
Kumar [281129] . 

In order to choose, for every substance, the most suitable equations of state 
to predict the minimal nucleation radius, we use the planar interface surface 
tension values in the following way: by following Cahn and Hilliard, we fix 
C (cf. the next subsection) while among the considered equations of state, 
we choose that which predicts values sufficiently accurately (for many 
substances this estimate holds within 90% of the experimental value). Finally 
we can predict the thickness of the interface and the surface mass density. We 
note that these data are relatively scarce especially the latter [T][rT] . 
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5. 1 The planar interface surface tension. 



Substance 


texp 


exp jdyn-i 
U °° I cm ' 


vdw \dyn-i 
00 L cm ' 




Ar 


~ 0.5 


13.1 


23.9 


3.24 


C0 2 


0.95 


1.16 


0.93 


1.51 


N 2 


~ 0.6 


9.80 


8.25 


1.50 


o 2 


~ 0.6 


13.2 


11.4 


1.46 


H 2 


~ 0.6 


2.09 


2.29 


1.34 


H 2 


0.5 


68.5 


63.8 


1.83 


NH* 


~ 0.7 


23.4 


26.3 


1.73 


CH3OH 


~ 0.6 


20.1 


20.1 


1.91 


C2H5OH* 


~ 0.5 


22.7 


24.7 


1.70 


c 2 h; 


~ 0.6 


16.2 


15.3 


1.50 


CeHia 


0.5 


25.5 


16.5 


1.63 


CCU 


~ 0.5 


29.8 


18.6 


1.64 


CF 4 


~ 0.9 


5.0 


0.87 
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Table 3. Planar interface surface tension for some substances evaluated in the 
case of van der Waals equation of state and the corresponding experimental value. 



c 



In the last column of the Table is reported the scaled parameter a exp — ~J^~^~- 

For experimental planar interface surface tension cf. |19|26j . The upper part of the 
Table shows substances for which Lennard-Jones radius r LJ is known with good 
approximation and so it is used to calculate doo by evaluating £ from (11/7) 5 r LJ . 
In the lower part of the Table (down the double line) the covolume radius r co is 
used to evaluate £ and then 0^. In the cases marked by * a better fitting was made 
by taking £ = 1.25 r co . 

In Table 3 we show collected experimental data for different substances. Values 
for the planar interface surface tension are shown together with the values 
evaluated by second gradient theory in the case of the van der Waals equation 
of state. The theoretical predictions were made by using in some cases the 
Lennard-Jones relation between r = r LJ and £ while in others we have simply 
put £ = r = r co , where r co is the covolume radius (we note that in general 
r LJ > r co but their ratio is lower than (11/7)2 ~ 1.25). 

The prediction of the second gradient theory is generally within 20% of the 
experimental value: this can be considered as a good result, comparable with 
the non-local model approach (cf. [8] and [IE]). In many cases our predictions 
are very accurate. We note however that are remarkable exceptions: 1) the Ar 
value is larger by about a factor two with respect to the experimental value; 
the same holds true in respect of the surface tension for other noble gases; 2) 
the values for C 6 H 12 and CC1 4 are smaller almost by a factor of two (about 
60% of experimental values), the same is true in minor form also for other 
carbon compounds such as C4H10 and C3H8 and most dramatically for CF4, 
where the ratio is about 1/6. 

The case of Argon and the other noble gases can be understood if we look to 
the last column of Table 3 in which we report the scaled experimental value of 

a exp = representing the normalized value of the van der Waals constant 

M. RT C 

a. For all substances a exp is between 1 and 2 whereas for Ar it is larger than 
2 (about 3.2); the same is true for other noble gases. This large value implies 
that the values of surface tension evaluated with the 'van der Waals theory' 
are generally far from the experimental values. This is a known result for Ar 
[T6] : the second gradient theory is insufficient to reproduce the experimental 
value for the surface tension, as higher-order gradients are necessary in the 
constitutive law for free energy to give a correct predictioi£0 (cf. [TB]). 

As for point 2) we could increase the values of the surface tension, by multiply- 



We note that values of a exp higher than for noble gases are obtained for water 
and methyl alcohol, this is reflected, (cf. Table 3, caption), in the use of r co rather 
than r LJ in the reported theoretical predictions. 
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ing r co by 1.25, i.e., the Lennard- Jones factor. This trick works effectively for 
ethyl alcohol, ethane and ammonia, but for the other substances the theoret- 
ical surface tension remains too small when compared with the experimental 
value (we note also the case of CO2 in which the theoretical value is 79% of 
the experimental value). So we conjecture that other equations of state, that 
generally have a larger value of o"oo (cf. Fig. 2), would describe carbon com- 
pounds better than the van der Waals equation. For the sake of brevity we do 
not present the theoretical data; we limit ourselves to saying that, by using a 
different equation of state, the values of can be predicted with the same 
precision as shown in Table 3 for substances that are well described by the 
van der Waals equation. 



5.2 The minimal nucleation radius and the other surface quantities 




Fig. 5. Minimal nucleation radii vs reduced temperature in units of C in the interval 
[0.5, 1.0[ for the following equations of state: van der Waals □; Rocard 3 O; Rocard 
4 *; Berthelot x; Peng-Robinson A; Eberhart-Schnyders +. State equations Rocard 
1 and 2 give minimal nucleation radii between Rocard 3 and van der Waals cases. 

In Figure 5 we display the normalized minimal nucleation radius in units of £ 
for the equations of state in Table 1 versus reduced temperature in the interval 
[0.5, 1.0 [. The interval chosen is the most relevant for many substances because 
most of the experimental data on surface tension exist in this range of reduced 
temperature (cf. [i~9T2"6] ). At the lower temperatures the minimal radius tends 
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to be a constant of approximate value C; the exact values could be retrieved 
in principle by the method described in Section 4, but the strong stiffness of 
the problem, induced by the rapid decrease to very small values of the gas 
phase density, limits the precision of the method. 



Substance 


Rm [A] 


Sm [A] 


p™ 10- 7 ■ [g/cm 2 ] 


CO2 


22.3 


9.52 


9.4 


N 2 


19.7 


8.47 


5.7 


O2 


18.4 


7.86 


7.4 


H 2 


15.6 


6.67 


0.4 


H 2 


15.5 


6.63 


4.6 


NH* 


20.8 


8.87 


4.5 


CH3OH 


20.2 


8.62 


5.0 


C2H5OH* 


27.3 


11.6 


6.9 


C 2 H* 


24.9 


10.6 


4.7 



Table 4. Minimal nucleation radii for bubbles of some substances evaluated at 
t = 0.85 in the case of van der Waals equation of state. In the other two columns 
are listed the predicted values of thickness and surface mass density at minimal 
nucleation radius at the same reduced temperature. The meaning of the division of 
the Table in upper and lower part and the * near some substance is the same than 
in Table 3. 

We see that the Berthelot equation of state and the modification proposed 
by Eberhart-Schnyders [JT] predict radii rather smaller than other equations 
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of state. The Peng-Robinson equation of state gives radii near to the van 
der Waals model at high values of t and near the Berthelot model values at 
about t = 0.5. The family of all equations of state proposed by Rocard give 
radii very near to the van der Waals value (for this reason only two of these 
data are reported in Figure 5, in order to leave the figure sufficiently clear). 
Obviously, when the critical temperature is approached the minimal radius 
and the thickness tend to oo. 



Substance 


State Equation 


Rm [A] 


Sra [A] 


p™ 10- 7 ■ [g/cm 2 ] 


CH 4 


Ebcrhart-Schnyders 


13.1 


5.70 


2.0 


C3H8 


Rocard 4 


22.1 


9.54 


4.6 


C 3 H 8 


Ebcrhart-Schnyders 


18.1 


7.86 


3.8 


C4H10 


Rocard 4 


25.1 


10.8 


5.4 


C4H10 


Eberhart-Schnyders 


20.5 


8.90 


4.5 


CeHia 


Rocard 4 


26.3 


11.3 


6.8 


CeHia 


Peng- Robinson 


23.6 


10.2 


6.2 


CCU 


Rocard 4 


26.1 


11.2 


13.8 


CCI4 


Peng- Robinson 


23.4 


10.1 


12.6 



Table 5. Minimal nucleation radii for bubbles of some substances (carbon com- 
pounds) evaluated at t = 0.85 in the case of some other equations of state. These 
will be chosen to reproduce the planar interface surface tension experimental values 
at least to 80%. In the other two columns are listed the predictions for thickness and 
surface mass density at minimal nucleation radius at the same reduced temperature. 

In Table 4, the minimal nucleation radii at t = 0.85 are listed for those sub- 
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stances whose planar interface surface tension is well described by the van der 
Waals equation of state. In the same table, we have added data on the thickness 
of the interface and the surface mass density for the substances cited therein. 
In Table 5, further theoretical predictions are given, always at t = 0.85, for 
substances whose planar interface surface tension is well described by other 
equations of state. Thus we show that radii, thickness and surface mass den- 
sities can be calculated very easily by using the results found in [13] . We note 
that the minimal radius ranges typically from 10 to 30 A so that it is just 
a few molecular radii. For example the molecules of CeH 12 have a radius of 
about 5 A; at t = 0.85 the ratio between this value and predicted minimal 
radius is about 4.7, while at t = 0.5 the minimal nucleation radius is ~ 9.8 
A and the ratio is about 2. 

Examples of predictions of the equilibrium radius as a function of R (or /i^) 
by means of classical theory are given by van Carey [UJ. However he limits 
himself to the case of 'large bubbles', the only ones for which classical theory is 
valid. The experiment of Fisher and Israelachvili [18] needs to be reinterpreted 
because the Laplace formula, on which is based the Kelvin relation used to 
extrapolate the radius, is not valid near the minimal nucleation radius |12j . 
Moreover, as noted above, their results can be compared with a process of 
homogeneous nucleation only after a thorough analysis of the effect of their 
experimental apparatus. We note that as cited above CeHi2 has a predicted 
minimal nucleation radius of about ~ 9.8 A at the temperature in which the 
experiment was performed (t ~ 0.5). It is larger than the radii reported by 
Fisher and Israelachvili (we also note that, on changing the equations of state 
or the evaluation of C, the minimal radius always remains greater than their 
data). On the other hand, data on homogeneous nucleation for very small radii 
exist for droplets (cf. [28i;29]) and comparison with our data shows that the 
nucleation radii of bubbles we predict are in some cases very close to nucleation 
radius for the droplets; for example, in the case of water, the Lamer and Pound 
data are in the range 8.0-8.9 A at t ~ 0.5, and our prediction is R m ~ 8.3 
A; the data for equilibrium radii reported by Kumar for the CCL at t ~ 0.5 
range between 12 — 17 A, these are very close to those we calculated which 
range (at t = 0.5) between 12 — 14 A, though this range of variation is due to 
the used equation of state rather than to the supersaturation ratio. 

Moreover we predict that the thickness of the interface is formed by few molec- 
ular layers [16] in the range from minimal nucleation radius to saturation. A 
new insight of the physical problem is obtained when observing that at small 
radii the interface constitutes a large part of the microscopic bubbles and 
cannot be confined in any Gibbs-Tolman idealized zero thickness layer. Typi- 
cally the numerically evaluated ratio between the equivalent thickness at the 
minimal radius and the minimal radius itself is of the order of 0.4 for all the 
equations of state studied here. Finally we note that the surface mass density 
values (see Tables 4 and 5) are of the order of magnitude predicted by Alts 
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and Hutter in their theory on water pQ. 



6 Conclusion 



The first three Sections of this paper give a complete account of the results 
available in the literature concerning the theoretical treatment of the phe- 
nomenon of homogeneous nucleation and the dependence of surface tension 
on curvature. 

However recently (cf. [12]) a rational definition for surface tension and radius 
of curvature for spherical interfaces has been proposed in terms of a solution 
of the density field, using a shell-like approach to the definition of interface 
quantities. Moreover, thanks to the second gradient theory, a deeper under- 
standing of the relationship between the Laplace formula and the capillary 
constant A has become possible. In particular our investigations stem from a 
theoretical treatment in which: 

1) an expression for surface tension is used which accounts for curvature ef- 
fects better than those found in the literature; 

2) the mechanical pressure is distinguished from the thermodynamic pressure. 

In this way we obtained, using numerical methods, predictions for: 

1) the dependence of surface tension, equilibrium radius, thickness of the 
spherical interface and surface mass density on chemical potential; 

2) the values of the minimal nucleation radius, i.e., the minimum equilibrium 
radius possible for a bubble. Being founded either on statistical mechanics or 
on second gradient theories, our predictions have a clear meaning in a range 
of temperature close to the critical point. 

Possible generalizations include the study of: 

(i) models in which A depends on p locally or globally; 

(ii) models in which higher density gradients appear, or a fully non-local ap- 
proach via integral equations; 

(iii) models suitable to describe the dynamics of phase transitions. 
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